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Proton-proton collisions at new high energies (^/s =2.36 and 7 TeV) at LHC resulted into greater mean multiplicities ((«)) of 
charged particles in the mid-rapidity region than estimated ones by different models and event generators. Another significant 
observation in multiplicity data is the change in slope in the distribution of primary charged hadrons in symmetric pseudorapidity 
interval \ri\ <2.4. The change is most prominent with data at yfs - 7 TeV. These new observations merit further studies. We 
consider a two-component model of particle production to analyze multiplicity distributions of charged hadrons from proton- 
proton collisions at centre-of-mass energies ^fs - 0.9, 2.36 and 7 TeV in symmetric pseudorapidity intervals \ri\ of increasing 
width around the centre-of-mass pseudorapidity //„„ - 0. The model, based on quantum statistical (QS) formalism, describes 
multiplicity distribution by convolution of a Negative Binomial Distribution (NBD), representing a chaotic component, and a 
Poisson Distribution (PD), representing a coherent component of particle productions. The behaviour of characteristic parameters 
of the model is followed by the LHC data, while a scaling law, involving information entropy in quantum statistical viewpoint 
and derived as a function of chaotic multiplicity obtained from the two-component model, is not obeyed by the data, satisfactorily. 
An attempt to match the measured multiplicity distributions and suggested convolutions with values of characteristic parameters 
extracted from the data confirms disagreement between the data and the model. 
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1. Introduction 

The probability distribution P„{n) of production of n parti- 
cles from collisions, has been the premier of basic observables, 
characterizing the final states of multi-particle production pro- 
cess in high energy physics experiments since the beginning of 
such studies. In recent times, the Large Hadron Collider (LHC) 
llll at CERN has taken high energy collisions in laboratories to a 
new energy domain, facilitating proton-proton (pp) collisions at 
centre-of-mass energies, yfs - 0.9, 2.36 and 7 TeV ||2> S B Si • 
The Compact Muon Solenoid (CMS) experiment at LHC has 
measured multiplicity distributions of primary charged hadrons 
for all the three LHC energies, available so far, in the mid- 
pseudorapidity (rj) region (where rj - -ln[tan{0/2)] and is the 
polar angle of the particle with respect to the counterclockwise 
beam direction) in five symmetric overlapping ?7-intervals \ri\ or 
rjc < 0.5, 1.0, 1.5, 2.0 and 2.4 |4] around the centre-of-mass 
pseudorapidity {rjcm = 0). The measured mean multiplicities 
at the new LHC energies ( =2.36 and 7 TeV) in the range 
of mid-pseudorapidity, have been found to be largely under- 
estimated by existing models (event generators like PYTHIA, 
PHOJET etc). Other experiments at LHC corroborate the find- 
ing, generating strong motivation to look into the LHC data 
from different approaches. Another significant observation by 
the CMS experiment is the change in slope in the distribution of 
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primary charged hadrons in pseudorapidity interval I?;! < 2.4 at 
^fs -1 TeV. Appearance of a change in slope in the measured 
distribution indicates to the possible existence of more-than- 
one process or source of particle production. In this article, we 
present results of our analysis of multiplicity distributions data 
from collisions at LHC, in the light of a two-component 
model 16t] formulated in a quantum statistical approach. We 
extend our study further to check the validity of a scaling law 
involving information entropy in quantum statistical point of 



2. Background 

In the last few decades, with the advent of accelerator tech- 
nology, collider facilities, capable of delivering higher and higher 
centre-of-mass energy ( ^fs), could be made available for proton- 
proton (pp) and proton- antiproton (pp) collisions. While the 
Intersecting Storage Ring (ISR), CERN facilitated pp collisions 
at range of 23.9 to 62.2 GeV |7], the Super Proton Syn- 
chrotron (SPS), CERN had pp colHsions at = 200 to = 
900 GeV [8]. In the Tevatron at Fermi lab, the energies of pp 
collisions were = 540 GeV to ^fs - \ .% TeV |9] and finally 
at the Large Hadron ColUder (LHC), CERN, energy of colli- 
sions of protons have reached yfs as high as 7 TeV flUlsl]. 
With all these data along with data at the fixed-target (non- 
collider) experiments at pre-ISR period, a rich set of experimen- 
tal data on multiplicity distributions from pp (pp) collisions is 
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now available for a wide range of a/s for a comprehensive and 
systematic study in different theoretical and phenomenological 
formalisms for better understanding of multiparticle production 
mechanism. The two-component model [6] of particle produc- 
tions is one such formalism which has been thoroughly used in 
analyzing experimental data of pp (pp) collisions at energies 
available up to the SPS. 

Of the statistical distribution functions, the Negative Bino- 
mial Distribution (NBD), 



P(n, <n>, k) 
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has been the most successful one in describing the probabil- 
ity distributions of final state charged particles from pp (pp) 
collisions in the discussed energy domain. At the lower part 
of the energy range, the multiplicity distributions of final state 
charged particles in the full phase space could be described by 
binomial distribution (when the parameter k in equation- 1 is 
negative and an integer). At around Vs = 5 GeV, the distribu- 
tion of produced charged particles turned broader and started 
following Poisson distribution (k is infinite). Above ^/s - 30 
GeV, the NBD matched fairly well with even broader distribu- 
tions in the full phase space (pseudorapidity space) data up to 
VJ = 540 GeV at SPS. At = 900 GeV at SPS, when a 
shoulder-like structure appeared in the tail of the multiplicity 
distributions, a single NBD failed in matching with the data at 
large pseudorapidity intervals. At 900 GeV and also at 1.8 TeV 
at Tevatron, the NBD function was successful only in the re- 
stricted mid-pseudorapidity (\ri\ < 0.5) region. The distribution 
in large 77-interval could be reasonably described with the sum 
of two NBDs. 

3. Motivation 

At the LHC energies, multiplicity distributions in non-single 
diffractive (NSD) inelastic proton-proton collisions have been 
measured and reported by different experiments, in different 
kinematic ranges, depending on capability of respective detec- 
tor setup in terms of geometrical acceptance, detection efficien- 
cies etc. A Large Ion Collider Experiment (ALICE) at LHC 
has measured primary charged particles at = 0.9 and 2.36 
TeV in the mid-// region in three overlapping //-intervals \r]\ < 
0.5, 1.0 and 1.3 [2]. At = 7 TeV, instead of NSD inelas- 
tic events, ALICE analyzed ^ an event class requiring at least 
one charged particle in I77I < 1 and measured multiplicity dis- 
tribution in that 77-interval only. The ATLAS experiment at 
LHC has measured charged particle multiplicities for different 
event classes characterized by different lower cuts on the num- 
ber of charged particles (nci, < 1,2 and 6) in diff'erent kinematic 
ranges (pr > 100 MeV, 500 MeV in |?7| <2.5). The distributions 
measured by ALICE in the three |?7|-intervals, I;;! < 0.5, 1.0 and 
1.3, at the two energies, = 0.9 TeV, V* = 2.36 TeV have 
been reported to match fairly well with NBD. The NBD fit to 
the distribution at ^fs -1 TeV, measured by ALICE, has been 
reported to be slightly underestimating the data at low multi- 
plicity (n < 5) and slightly overestimating the data at high mul- 
tiphcity (« > 55). Measurements of multiplicity distributions 



in a wider phase space and in all the three energies by CMS 
experiment, reveal interesting features. Presenting |4] multi- 
plicity distributions, without fitting to any distribution function, 
CMS reported a change of slope in P„ for n > 20 in its largest 
77-interval of I77I < 2.4. This feature becomes more pronounced 
with increasing yfs. 

The observed change in slope in the distributions as mea- 
sured by CMS could be attributed to the existence of more than 
one kind of source or process of particle production. There have 
been several proposals of models |0, T^, U, 12, 13] involving 



multi-particle production from more than one process / source. 
In Ref.|10], it is the weighted superposition of two multiplic- 
ity distributions (each assumed to follow NBD); one due to soft 
events (without mini-jets) and the other representing semi-hard 
events (with mini-jets). In the multiparton interactions model 
described in Refs. ll iL 11211 . the soft component, which consti- 
tute the bulk of events in the final state, corresponds to sin- 
gle parton-parton collisions in the framework of dual parton 
model and produces KNO [14,] distribution. The other process, 
which seems to be superimposed on top of the KNO producing 
process, involves two or more independent parton-parton colli- 
sions. Ref. lil3i] is a two-component dual parton model (DPM), 
the soft process is described by the supercritical Pomeron and 
the hard component is described by perturbative QCD. This 
two-component DPM includes diffl^active processes also. In 
spite of all these efforts, the change in slope in distributions 
is not yet a fully understood phenomenon. 

We recollect that the distribution at ^fs = 900 GeV in broad 
phase space has been explained by sum of two NBDs. Also, 
there has been agreement (l^ on the soft parts (as measured in 
UA5 and E735 experiments) obtained by the models described 
in Ref.|[l3] and Ref.liij]. But, at this point, it is important to 
consider relative position of the change in slope in the distribu- 
tions of data at -v^ = 900 GeV at SPS and a/s = 7 TeV at LHC. 
The change in slope appeared in SPS data at a higher multiplic- 
ity (n) and a sum of two NBDs (representing two broad distribu- 
tions) could describe the data satisfactorily. For the LHC data, 
the most prominent change in slope appears in the multiplicity 
distribution in the symmetric pseudorapidity interval, 77c < 2.4 
at ^/s -1 TeV at much lower multiplicity (77)20), indicating that 
the measured distribution could be due to contributions from a 
broad distribution and a narrow distribution. 

Such a scenario is provided by a two-component model 

in quantum statistical approach, where a weighted 
average of a NBD (representing a broad distribution) with k - \ 
due to a chaotic source and a Poisson distribution (representing 
a narrow distribution) due to a coherent source results into the 
final distribution. We follow this two-component model in QS 
viewpoint to analyze the LHC data. 



4. Methodology 

In the two-component model 
mean multiplicity is given by: 

<7?) = {rich) + (rico) 



in QS approach, the total 
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where {rico) is the coherent component and {rid,) is the chaotic 
component of the mean of total multiplicity. 

In quantum statistics, a completely chaotic source produces 
particles following a negative binomial distribution given by 
equation- 1 (n in equation- 1 is to be read as rich in the present 
formalism), where k is the number of cells in the phase space or 
the number of independent quantum states. On the other hand, a 
completely coherent source gives rise to a Poisson distribution: 



P{{nco),kp)^ 



(3) 



The mean of the chaotic component of total multiplicity, 
{rich) is obtained from the measured mean of total multiplicity 
and the second moment of the distribution by: 



{rich) = p{n) 



where 



p = [K{C2-il + ll{n))]] 



1/2 



(4) 
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is the chaoticity parameter, C2 = {rP')l{rif' is the second mo- 
ment of multiplicity distribution and K can have value either 
1 or 2. However, in narrow symmetrical rj intervals in data, 
only ^"=1 leads to physical solutions, in most of the cases. 
The chaoticity, p {0 < p < I) indicates measure of fraction 
of chaoticity involved in the source. For a completely chaotic 
source, p - I and for a completely coherent source, p - Q. 
Any value, in between, represents a convoluted distribution of 
the two. 

For our present study, we use the published |4] multiplic- 
ity distribution data from CMS experiments in five overlap- 
ping symmetric pseudorapidity intervals I;;! around the centre- 
of-mass pseudorapidity with width of the interval extended up 
to I77I < 2.4 for the three LHC energies, V* = 0.9, 2.36 and 7 
TeV. In some cases, we use published ALICE data in \ri\ < 
0.5, 1.0 and 1.3 for ^fs = 0.9 and 2.36 TeV for which the mul- 
tiplicity distributions in non-single diiTractive (NSD) inelastic 
proton-proton collisions have been measured and published by 
the experiment. 

5. Results and Discussions 

5.1. Behaviour of characteristic parameters of the model 

Following Fowler et al.|'6'], using equations-5 and 4 given 
in section-4 above, we calculate chaoticity {p) and mean of the 
chaotic multiplicity {{rich)) for all the multiplicity distributions 
in NSD pp collisions as reported|2, 4] by the two LHC experi- 
ments, CMS and ALICE. We plot the parameters in Fig.[T]and 
Fig. |2] respectively. In these plots, we include similar data for 

= 540 GeV of SPS also for comparison. The chaoticity de- 
picts fractional contribution of chaotic source in the total multi- 
plicity. As could be seen in Fig.[T] the chaoticity decreases with 
increasing width of the 77-interval, 77c for a given ^fs. Also, for 
a given width of rjc, the chaoticity increases with ^[s. 
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Figure 1: Chaoticity for = 0.54, 0.9, 2.36 and 7.0 TeV for different r]- 
intervals. Tlie error-bars, include both the statistical and the systematic uncer- 
tainties. Some of the error-bars are drawn inclined for clarity in presentation. 



It is worth mentioning at this point that, for the LHC data, the 
statistical and the systematic uncertainties as quoted with prob- 
ability distribution data |4fl have been added in quadrature 
and propagated all through the analysis. The error-bars associ- 
ated with the LHC-data points in the plots in all the figures of 
this article, therefore, include both the statistical and the sys- 
tematic uncertainties. The plots of chaotic multiplicity in Fig.|2] 
shows rapid increase in chaotic multiplicity with both the rjc and 

We also calculate the coherent component ({rico)) of mean 
multiplicity using equation-2 in section-4 and study the depen- 
dence of {rico) on rjc and ^fs. Our analysis reveals (as has been 
shown in Fig.O almost no change in {rico) with increasing ^fs 
in a given rjc and a very low rate of increase in (as com- 
pared to the rate of increase in {rich)) with increase in rjc for a 
given Vi, in the considered mid-pseudorapidity region. 
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Figure 2: Mean of chaotic multiplicity for = 0.54, 0.9, 2.36 and 7.0 TeV for 
different (^-intervals. Lines drawn by joining data-points to guide the eye. The 
error-bars include both the statistical and the systematic uncertainties. 
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Figure 3: Mean of coherent multiplicity for = 0.54 TeV, 0.9 TeV, 2.36 TeV 
and 7.0 TeV for different ;/-intervals. The en'or-bars include both the statistical 
and the systematic uncertainties. 



Figure 4: The information entropy, calculated with parameters of the two- 
component model in QS approach, as extracted from the LHC data do not ap- 
pear to obey the scaling law ll^ . 



Comparing observed dependences of the parameters of the model 
for the LHC data on rjc and -^/s, the particle production scenario 
at the considered energy range by the present approach could be 
summerized as follows: 

a) The chaotic source populates mainly the mid-pseudorapidity 
region and the fractional contribution of the chaotic source as 
compared to the coherent source decreases as one goes away 
from the centre-of-mass pseudorapidity {t]c„, - 0). 

b) The span of pseudorapidity of dominance of chaotic source 
increases with increase in ^fs. 

c) The increase in the total mean multiplicity («) with Vs in 
the mid-pseudorapidity region, in the considered energy range, 
is dominantly due to the increase in chaotic multiplicity {rich} 
with yfs. 

In general, the behaviour of the parameters extracted from 
the LHC data match with characteristic features of the model 
in a similar way as has been observed in case of SPS energy 
data. It is naturally logical to analyze the LHC data in terms 
of a scaling law U3], involving chaotic multiplicity, which has 
been reported to hold for pp or pp data for a wide range of 
energies at ISR and SPS. 



It may be noted that the discussed scaling law is not a charac- 
teristic feature of the two-component model. The scaling law 
was developed for information entropy in quantum statistical 
point of view, modifying the entropy scaling 1 19] found to obey 
experimentally measured multiplicity distribution in hadronic 
interactions for wide range of ^/s, including ISR and SPS data. 
The response of the LHC-data to the scaling law is worth ob- 
serving. We study the scaling law with the CMS data of pp 
collisions at = 0-9, 2.36 and 7 TeV along with the SPS 
data of pp at ^fs -540 GeV. The data points along with respec- 
tive experimental uncertainties are plotted and are joined with 
straight lines for a given yfs to guide the eye. As it is clear from 
the Fig.m the data points of different yfs, do not follow a com- 
mon single curve. Moreover, data-points of different energies 
follow distinctly separate lines. This observation indicates to 
the deviation of the entropy-scaling at LHC energies. 



5.2. Scaling of information entropy 

We analyze the LHC data of NSD pp collisions up to ^fs 
= 7 TeV, as measured by the CMS experiment to test the valid- 
ity of a proposed scaling law namely the scaling of information 
entropy [18] involving application of the discussed two compo- 
nent model. The information entropy is a function of the chaotic 
multiplicity (< n^, >), which is a function of symmetric pseu- 
dorapidity interval, rjc and the centre-of-mass energy, ^/s. The 
entropy for chaotic multiplicities in symmetric pseudorapidity 
intervals 77c is given by: 

S(t]c, ^/s) = (< «cA > +!)/«(< nch > +1) 

- <nch> In < nch > (6) 

In Fig. m we plot S lri„,ax as a function of ^ = rjc/rimax where 
ill, 

Imax = ln[( yfs - 2mn)/m„] (7) 



5.3. Test of agreement between the model and the data. 

At this stage, one might be interested to know how well 
the discussed two-component model describes the measured 
multiplicity distributions. To check the agreement between the 
model and the data, we chose multiplicity distributions in the 
pseudorapidity interval, I77I <2.4 for the fact that the change in 
slopes in distributions have been observed by the CMS collabo- 
ration for in the same 7;-interval with most pronounced structure 
appearing at yfs - 7 TeV and that the motivation of this work 
has been the understanding of change in slope of the multiplic- 
ity distribution by the two-component model. 

We fit the multiplicity distributions of primary charged hadrons, 
as measured by the CMS experiment, at energies ^/s - 0.9, 2.36 
and 7 TeV in the pseudorapidity interval {rjc <2.4) around the 
mid-?7 with NBD. For clarity in presentation, plots of three en- 
ergies are given separately in Fig. |5] |6] |2l The lower panel of 
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Figure 5: Primary cliarged hadron multiplicity distributions for <2.4 for ^fs 
= 7.0 TeV. The continuous line corresponds to fit to a NBD function and the 
dotted line corresponds to the convolution of a NBD and a PD with parameters 
derived from the two component model as described in Ref.[6]. The lower 
part of the panel contains ;t'^/dof as well as the plots of residual analysis to test 
goodness of the fits. The error-bars shown with the data-points include both the 
statistical and the systematic uncertainties. 



Figure 6: Same as Figure - 4 for ■\/s = 2.36 TeV. Lines of dift'erent styles 
correspond to different fits similar to those as described in the caption of Figure 
-4. The lower part of the panel contains x^/^of ^ well as the plots of residual 
analysis to test goodness of the fits. The error-bars shown with the data-points 
include both the statistical and the systematic uncertainties. 
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each of these figures contains information on goodness of re- 
spective fits. Interestingly, all three distributions, including the 
one for rjc <2.4 at = 7 TeV, where a prominent change in 
slope appears, match well with NBD as it is evident from the 
values of the x^/dof, quoted in respective figures. In the same 
figures the dotted lines correspond to the fit with the convolu- 
tion of a NBD (with k=l) and a PD with parameters extracted 
from the data by the discussed two component model. 
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model of particle production based on quantum statistical for- 
malism. Though the behaviour of characteristic parameters of 
the formalism is found to be consistent with the LHC data, a 
scaling law of quantum statistical information entropy, calcu- 
lated from chaotic multiplicity, extracted from the LHC data by 
the model, is violated. 

A test of agreement between the model and the data by 
fitting the multiplicity distributions with convolutions of func- 
tions as suggested by the model and using values of the char- 
acteristic parameters evaluated from the data reveals that the 
model indeed fails to describe the data satisfactorily. On the 
other hand, the commonly used single NBD function fits better 
than the convolution suggested by the model. Our analysis thus 
shows disagreement between the LHC data and the discussed 
two-component quantum statistical model in its present form. 
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Figure 7: Same as Figure - 4 for ^fs = 0.9 TeV. Lines of different styles cor- 
respond to different fits similar to those as described in the caption of Figure 
-4. The lower part of the panel contains ;(f-/dof as well as the plots of residual 
analysis to test goodness of the fits. The error-bars shown with the data-points 
include both the statistical and the systematic uncertainties. 



As can be seen by visual inspection of the plots (with zoomed 
insets) as well as from the corresponding plots of residual anal- 
ysis and the values of x^l^o^, the multiplicity distributions fit 
better with single NBD than with the convolution of a NBD 
and a PD in terms of the discussed two component model, in- 
dicating poor agreement between the model and the LHC data 
of multiplicity distribution. The residual is defined as the 
difference between a data-point and corresponding fit-value. 



6. Summary and Remarks 

We have studied multiplicity distributions of primary charged 
hadrons produced in proton-proton collisions of non-single diffrac- 
tive class of events at LHC energies, in terms of a two-component 
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